Search for inelastic dark matter with the CDMS II experiment 
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Results are presented from a reanalysis of the entire five-tower data set acquired with the Cryo- 
genic Dark Matter Search (CDMS II) experiment at the Soudan Underground Laboratory, with 
an exposure of 969kg-days. The analysis window was extended to a recoil energy of 150 keV, and 
an improved surface-event background-rejection cut was defined to increase the sensitivity of the 
experiment to the inelastic dark matter (iDM) model. Three dark matter candidates were found 
between 25keV and 150 keV. The probability to observe three or more background events in this 
energy range is 11%. Because of the occurrence of these events the constraints on the iDM param- 
eter space are slightly less stringent than those from our previous analysis, which used an energy 
window of 10-100 keV. 



PACS numbers: 95.35.-|-d, 29.40.Wk, 95.30.-k, 95.30.Cq 

Cosmological observations [l| strongly suggest that 
nonluminous, nonbaryonic matter constitutes most of the 
matter in the Universe. This dark matter should be dis- 
tributed in dark halos of galaxies such as the Milky Way, 
enabling the direct detection of the dark matter parti- 
cles via their interactions in terrestrial detectors [2] . The 
movement of the Earth around the Sun would provide 
an annual modulation of the counting rate, caused by 
the change in the relative velocity of the dark matter 
particles with respect to the earthbound target 

The DAMA collaboration claims the observation of 
such a modulation in two different Nal(Tl) scintillation 
detector arrays 0, [El- The observed signal is in the 



2-6 keV electron-equivalent energy range with a period- 
icity of 0.999±0.002 years and a phase of 146±7 days [6]. 
The observed modulation signature is consistent with the 
expected signature of galactic dark matter particles in- 
teracting in a terrestrial detector. Other experimental 



results [7H12||. however, are inconsistent with the inter- 
pretation of the DAMA result as a si gna l from weakly 
interacting massive particles (WIMPs) l3- 13] elastically 
scattering off nuclei. 

Inelastic dark matter (iDM) scattering has been pro- 
posed as a way to resolve this tension [17j . The inelastic 
scenario assumes that WIMPs (x) can only scatter off 
baryonic matter {N) by transition into an excited state at 
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a certain energy above the ground state (x iV — >• x* -/V), 
while elastic scattering is forbidden or highly suppressed. 
There is a minimal velocity required to produce recoil 
energy Eji in such an inelastic scatter, 

1 / tunEr \ 

Vmin = ^==== +<5 , (1) 

where m^v is the mass of the target nucleus, /i is the re- 
duced mass of the WIMP-nucleus system, and 6 is the 
WIMP-mass splitting; ^ = keV is equivalent to elas- 
tic scattering. If Er is too small or too large, Wmin is 
above the cutoff imposed by the galactic escape velocity, 
and the event cannot occur. Important consequences of 
this model for direct detection experiments are differen- 
tial rates that peak at tens of keV recoil energy, and a 
significant suppression of the recoil spectrum at low recoil 
energies. In addition, the annual modulation signature 
is significantly enhanced, because of the increased de- 
pendence on the high- velocity tail of the WIMP-velocity 
distribution, which in turn is due to the larger minimal 
velocity (see Eq. ([T])). Therefore, the iDM scenario is also 
particularly sensitive to the escape-velocity cutoff in the 
WIMP-velocity distribution. Finally, it is important to 
note that the scattering rate is enhanced for heavy target 
nuclei (e.g. Xe and I). 

The Cryogenic Dark Matter Search (CDMS II) experi- 
ment operated in the Soudan Underground Laboratory a 
total of 19 Ge (-230g each) and 11 Si (~105g each) de- 
tectors at a temperature of ^^40 mK [3) [3| • These semi- 
conductors were stacked into five towers (T1-T5) with six 
detectors (Z1-Z6) each. They were instrumented with 
four channels of superconducting transition-edge sensors 
on the top side to detect phonons and two concentric 
electrodes on the bottom side to detect ionization. The 
primary ionization signal was read out by an inner elec- 
trode covering ^^85% of the detector surface. A thin 
outer electrode served as a guard ring to identify and 
reject events at the edge of the detector, which was sub- 
ject to higher background and reduced charge collection. 
The recoil energy was reconstructed from the phonon and 
the ionization signal [20] . The ratio of ionization to re- 
coil energy ( "ionization yield" ) was lower for nuclear re- 
coils, produced by WIMP candidates, than for electron 
recoils, caused mostly by background photons. Fewer 
than 10""' of the electron recoils in the bulk of the de- 
tector were misidentified as nuclear recoils. The main 
source of misidentified electron recoils were events with 
interactions in the first few /itm of the detector surfaces. 
Because of incomplete charge collection these events had 
reduced ionization yield, and occasionally the reduction 
was severe enough to mimic a WIMP-nucleus interac- 
tion. The phonon signals of these surface electron-recoil 
events had faster-rising pulses than bulk nuclear recoils 
and occurred closer in time to the more prompt ioniza- 
tion pulses. As discussed in detail in this paper a cut 
based on these timing parameters was employed to re- 



ject interactions at the detectors' surfaces. Misidentified 
surface events constituted the dominant background for 
the CDMS II experiment, while the neutron background 
from cosmogenics and radioactive processes was much 
less significant. 

Initial constraints from CDMS on the iDM model in- 
terpretation of the DAMA claim were set using a recoil- 
energy range of 10-100 keV 0. This paper presents a 
dedicated iDM analysis of the entire CDMS II five-tower 
data set, taken during two periods of stable operation be- 
tween October 2006 and July 2007 (internally denoted as 
runs 123-124) [2lj, and four periods between July 2007 
and September 2008 (internally denoted as runs 125-128) 
0. Note that the constraints on the WIMP-parameter 
space shown in Q were a combination of the final results 
from all data sets taken at the Soudan Underground Lab- 
oratory, which, however, were analyzed separately. In 
particular, the surface-event rejection cuts, as discussed 
below, were set at fixed backgrounds for runs 123-124 
and runs 125-128 separately. For the analysis presented 
here, the whole acquired data were combined in advance 
and surface-event rejection was based on the entire data 
set. There were two main reasons for performing this re- 
analysis. The iDM parameter space allowed by our pre- 
vious analysis (see Fig. 4 of (ll) includes WIMP masses 
mw ~ lOOGeV/c^ and mass splittings 6 ~ 120 keV. As 
shown in Fig. [l] these parameters result in a significant 
expected rate above our previous analysis upper limit of 
100 keV, so a simple extension to 150 keV increases the 
expected sensitivity. Moreover, the expected rate drops 
to zero for low recoil energies, in contrast to the elastic- 
scattering case, obviating the need for a low threshold. 




Recoil energy [keV] 

FIG. 1. (color online). Differential recoil spectra in a Ge 
target for a WIMP-mass splitting S of 120 keV and a few rep- 
resentative WIMP masses mw- For comparison tiie spectrum 
for a WIMP with a mass of 120 GeV/c^ assuming elastic scat- 
tering {S — keV) is also shown (black/solid). The spectra 
are normalized to unity in the 10-150 keV recoil-energy range. 
The vertical lines denote the analysis threshold at 10 keV, the 
lower boundary for the setting of the surface-event rejection 
cut at 25 keV, and the upper analysis limit from our previous 
analysis at 100 keV [3| . See text for details. 
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Since most of the dominant surface-event background oc- 
curred at energies just above our lOkeV threshold plj . 
where no iDM signal is expected, the sensitivity could 
be further improved by redefining a looser surface-event 
rejection cut based upon the estimated background with 
recoil energy between 25keV and 150 keV, while leaving 
the lower boundary for the analysis at lOkeV. Thus, a 
significant number of surface-background events was ex- 
pected in the 10-25 keV range, which, however, had only 
a minor effect on the results in the parameter-space re- 
gion of interest {mw ~ 100GeV/c^ S ~ 120 keV). 

The same data-quality selection cuts used in previous 
analyses 0, [HI for ensuring detector stability and re- 
moving periods of poor detector performance, e.g. due 
to insufficient neutralization, causing incomplete charge 
collection owing to impurities in the detector crystal, re- 
sulted in a total Ge exposure of 969 kg-days for this re- 
analysis. The Si detectors were omitted due to their 
lower sensitivity to inelastic scattering. Because both 
data sets had already been analyzed, this analysis was 
not "blind" . However, the analysis was performed in a 
similar manner to minimize bias: selection criteria and 
background estimates were defined and evaluated using 
only WIMP-search data outside the signal region and 
calibration data. 

The detectors were exposed to gamma rays from -'^■^'^Ba 
and neutrons from ^^^Cf at regular intervals to calibrate 
their response and define criteria for data-quality cuts 
and the WIMP-acceptance region. The latter was de- 
fined to be the zL2a band around the mean nuclear-recoil 
ionization yield in the yield versus recoil-energy plane. 
An illustration is given in Fig. [2j which shows ^^^Cf cali- 
bration data from a representative detector in one of the 
six data runs. 

In addition to the quality cuts, most of the selec- 
tion criteria for WIMP-nucleon interactions remained un- 
changed from the previous analyses 0,121}. This included 
the single-scatter cut, requiring there to be no signal ex- 
ceeding the phonon-noise level by more than Aa in any 
of the other 29 detectors; the ionization-based fiducial- 
volume cut, rejecting events near the edges of the detec- 
tors; and the muon-veto cut, demanding negligible co- 
incident energy deposited in the active muon veto sur- 
rounding the apparatus. 

Extending the analysis window to 150 keV was hin- 
dered by the fact that statistics from the ^^^Cf neutron 
source were low above ~100keV which can be seen in 
Fig. [21 Thus, we extrapolated the nuclear-recoil bands 
at higher energies from the fits below 100 keV. The ex- 
trapolation showed good agreement with Lindhard the- 



ory 
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23l | when statistics from all six runs were com- 
bined for each detector, and both the band locations and 
the nuclear-recoil cut efficiencies had only a minor energy 
dependence above ~25 keV. 

The surface-event rejection was based upon a "timing 
parameter" consisting of the sum of the rise time of the 
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FIG. 2. (color online). Ionization yield versus recoil energy 
of ^^^Cf calibration data from a representative detector in 
one of the six data runs. The black/solid lines represent the 
electron-recoil band around a yield of one and the nuclear- 
recoil band around 0.3. The black/dashed line denotes the 
mean of the latter band, while the similar but blue/solid line 
is the corresponding prediction from Lindhard theory [2^. [2^. 

largest phonon pulse and its delay relative to the ioniza- 
tion pulse. This timing cut was set in the 25-150 keV 
energy range using Ba and Cf calibration data. Since 
surface events in WIMP-search data did not have the 
same recoil-energy and ionization-yield distributions as 
in Ba calibration data 0], this cut was not expected to 
be optimal, although corrections based on WIMP-search 
multiple scatters were applied to the distributions to di- 
minish these differences. Thus, the cut performance had 
to be tested on WIMP-search data before "unblinding" . 
The cut setting and testing are discussed in more detail 
in the following two sections. 

Each detector had its own timing-parameter cut. We 
tuned the set of cuts to yield a given expected "leak- 
age" (number of background events) for the whole data 
set, while maximizing the signal for a WIMP of mass 
lOOGeV/c^ and a mass splitting of 120 keV. For each 
given expected leakage, using values in steps of 0.1 be- 
tween 0.1 and 1.5, we ran Monte Carlo simulations to 
find the average upper limit we could obtain if there 
were no true WIMP signal. For each expected leak- 
age, 10^ surface-event mock data sets were generated, 
each with number and energies of background events 
chosen randomly according to the given expected leak- 
age and the expected energy distribution as estimated 
from WIMP-search multiple scatters. As was to be 
done with the actual data, a 90% C.L. upper limit on 
the spin-independent WIMP-nucleon cross section asi 
was calculated for each mock data set, using the opti- 
mum interval method 2J| with the WIMP recoil-energy 
distribution [2^, [2^ for the chosen WIMP parameters 
{mw = 100GeV/c2, S = 120 keV). Figure [3] shows the 
mean upper limit obtained as a function of the expected 
leakage used in selecting the set of timing-parameter cuts. 
The timing-parameter cuts were finalized at the values 
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obtained for a fixed expected leakage of 0.6 events, where 
the minimum was reached. 

As explained above, the leakage value chosen for op- 
timizing cuts was not a sufficiently accurate estimate 
of the expected background. Thus, as with our earlier 
analysis [7|, we used WIMP-search data to improve our 
estimate of the expected leakage. We estimated the leak- 
age by multiplying the number of WIMP-search nuclear- 
recoil single scatters failing the timing cut by pass-fail 
ratios deduced from event samples which were assumed 
to resemble the population of background events. For de- 
tectors that were not located at the top or bottom of their 
towers (interior detectors), two classes of multiple-scatter 
events in the WIMP-search data were used independently 
to estimate the ratios, and therefore the expected back- 
ground: events with ionization yield within the nuclear- 
recoil band, and events in which a detector had yield 
just above or below the nuclear-recoil band (wide-band 
events). The latter class was defined to include events 
outside the ±2cr nuclear-recoil band that had an ioniza- 
tion yield above 0.1 and below the minimum of 0.7 and 
the lower boundary of the ±5(T electron-recoil band at 
the events' recoil energies. We also included two detec- 
tors at the bottom of their towers (end cap detectors) 
in this analysis. In this case, we treated surface events 
on the top (internal) and bottom (external) sides of the 
detectors separately. The pass-fail ratios of the internal 
sides were estimated from multiple-scatter events with 
ionization yield within the nuclear-recoil band, and those 
of the external sides, where tagging of multiple scatters 
was not possible, from single scatters within the wide- 
band region. In both cases, interior and end cap de- 
tectors, we applied appropriate correction factors to the 
pass-fail ratios from wide-band events to account for dif- 
ferences in timing performance between surface events 
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FIG. 3. Expected sensitivity of this analysis for a WIMP of 
mass lOOGeV/c^ and mass splitting 120 keV for different pre- 
defined surface-event leakages at which the timing cut could 
be fixed. The cut corresponding to the minimum of 0.6 events 
was chosen as the final cut for this analysis before looking at 
the WIMP-search signal region. 



within and outside the nuclear-recoil band. For the end 
cap detectors, additional correction factors were intro- 
duced to correct for differences in the single-scatter event 
fractions on the top and bottom sides. Systematic errors 
from the estimates of these correction factors, as well 
as from systematic differences in timing-cut performance 
between single and multiple scatters, were included in 
the leakage calculation [26| . Because of the low number 
of events passing the timing cut a dedicated Bayesian 
surface-event leakage estimate was applied [2^, adding 
another systematic error from the choice of prior distri- 
bution. The final background distribution obtained by 
combining the two estimates from the interior detectors 
with the estimate from the end cap detectors is shown in 
Fig.m It contains all statistical and systematic errors. It 
has a maximum around 0.6 events where the leakage had 
been fixed for the setting of the cut, while the median, 
which we use as the final background estimate, is slightly 
higher but agrees with this value within error bars: 

M25-i50kcV - 0.8t°j(stat)t°;3(syst) . (2) 

As expected, a similar estimate in the low-energy range 
from 10-25 keV resulted in a much higher number of ex- 
pected leakage events: 

Mio-25kcV = 5.7tl;l{sta.t)to^{syst) . (3) 

Figure [5] compares the final efficiency from this analy- 
sis, based on runs 123-128, with the efficiency from the 
previous analysis of runs 125-128 [?[• In both analyses 
the surface-event rejection cuts had roughly the same ex- 
pected leakage in the energy range the cut was defined on 
(10-100 keV for the previous analysis and 25-150 keV for 
the analysis presented here). Even though an exposure 
which was larger by a factor of 1.6 was considered for the 
setting of the timing cut, the final efficiency increased by 
a factor of ~1.5. This improvement in efficiency was pos- 
sible because we neglected background at energies below 
where a signal is expected from iDM. 

1.2| 1 1 1 1 1 1 1 1 
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FIG. 4. Distribution of the surface-event leakage estimate in 
the 25-150 keV energy range. See text for details. 
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FIG. 5. (color online). Comparison of the total nuclear-recoil 
efficiency obtained in this analysis (red/dashed) and from our 
previous analysis based only on runs 125-128 (black/solid) 01 . 
The latter is only defined up to 100 keV. Redefining the timing 
cut achieved an increase in efficiency by a factor of ~ 1.5. As 
in Fig. [1] the vertical lines denote the analysis threshold at 
10 keV, the lower boundary for the setting of the surface-event 
rejection cut at 25keV, and the upper analysis limit from our 
previous analysis at 100 keV. 

Neutrons, induced by muons and produced by radioac- 
tive processes within the experimental apparatus, consti- 
tuted an additional background vi^hich was indistinguish- 
able from a WIMP interaction in the detectors. Exten- 
sive simulations carried out with GEANT4 
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and 

FLUKA [29, 30] indicated that the neutron background 
in the 25-150 keV energy range inducted by muons is ex- 
pected to be 0.04to'o3(stat.), and the background from 
radioactive processes is estimated to be between 0.03 and 
0.06. The background between lOkeV and 25keV is pre- 
dicted to be 0.06l^;J5!^(stat.) from muon-induced neutrons 
and between 0.04 and 0.08 from radiogenic neutrons. 
These estimates are valid for the reanalyzed exposure 
and include cut acceptances. 

After "unblinding" , eleven events were observed within 
the acceptance region passing the surface-event rejec- 
tion cut, three within the 25-150 keV range and eight 
between lOkeV and 25keV. Figure 6 shows these can- 
didates, along with all other WIMP-search events in or 
close to the signal region, which pass all constraints ex- 
cept for cuts on the ionization yield and timing param- 
eter. As can be seen in Table HI the candidates are well 
distributed over the whole data-taking period and are 
spread over various detectors; though, six of the eleven 
candidates occurred in the two end cap detectors (T3Z6 
and T4Z6), where there was less shielding from back- 
ground, and where there was no detector below it to help 
reject background by detecting multiple scatters. It was 
verified that the performance of the experiment was sta- 
ble at the times during which the events occurred. The 
detectors in which the three candidates above 25 keV 
occurred are examined in more detail in Fig. [71 where 
normalized ionization yield, defined as the distance from 
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FIG. 6. (color online). Ionization yield versus recoil en- 
ergy for all low-yield WIMP-search single scatters from the 
combined five-tower data set (runs 123-128) passing all con- 
straints except for cuts on the ionization yield and timing pa- 
rameter. Most events fail the timing cut (green dots), while 
most of the few dozen that pass the timing cut ( x ) fall out- 
side the nuclear-recoil band (blue/solid lines). Eleven events 
pass all the selection criteria, with three occurring within 
the 25-150 keV range upon which the surface-event rejection 
cut was defined, and eight between the lOkeV recoil-energy 
threshold and 25 keV. The ionization-energy threshold is also 
shown (blue/dashed-dotted). This threshold and the shown 
nuclear-recoil band represent the exposure-weighted means 
over all runs and detectors. 

the nuclear-recoil band mean measured in units of stan- 
dard deviations given by the width of the band, is plotted 
against the timing parameter relative to the actual cut 
position. The black/solid line denotes the timing-cut 
boundary on the given detector, and the shaded box in- 
dicates the acceptance region. The top plot is for T4Z6, 
with a WIMP candidate at 37.3 keV and three additional 
candidates below 25 keV. T4Z6 was a detector at the bot- 
tom of its tower with reduced background rejection ca- 
pability. The middle plot shows events in T4Z2, where 
an event occurred just above the timing-cut boundary 
with a recoil energy of 73.3 keV. Finally, we show events 
from T1Z2 in the bottom plot with a candidate above the 
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TABLE I. Distribution of the eleven event candidates over 
detectors and time. The two marked events (*) are the candi- 
dates already found in our previous analysis of runs 125-128 
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FIG. 7. (color online). Number of standard deviations each 
event is away from the mean of the nuclear-recoil band (nor- 
malized ionization yield) versus timing parameter relative to 
the timing-cut position (shifted timing parameter) for the 
three detectors with WIMP candidates above 25keV. The de- 
tectors are (from top to bottom) T4Z6, T4Z2 and T1Z2 with 
candidate events at 37.3 keV, 73.3 keV and 129.5 keV, respec- 
tively. In detector T4Z6 three additional candidates occurred 
in the 10-25 keV range. The acceptance regions are indicated 
by the shaded boxes. 

analysis range from previous analyses at 129.5 keV. This 
event is far above the timing-cut boundary and would be 
rejected neither by the surface-event cut from the previ- 
ous analysis [vj, nor by a tighter timing cut tuned to a 
leakage as low as 0.1 (instead of 0.6) events, which was 



the most stringent timing cut we tested. No additional 
events appear in the signal region above 25 keV until the 
timing cut is loosened to an estimated surface-event leak- 
age of more than 2.0 events. 

The probability to observe three or more surface- 
leakage events between 25keV and 150 keV given the 
background distribution f{^) shown in Fig. 2] was cal- 
culated as 



fc=3 



e '^^'^ 
fc! 



(4) 



and yields 9%. Inclusion of the estimated neutron back- 
ground increases this probability to 11%, which is low but 
not negligible. Thus, this analysis does not constitute a 
significant detection of WIMP scattering. The probabil- 
ity for eight or more surface-background events between 
the lOkeV threshold and 25keV was calculated based on 
a background distribution obtained analogously to the 
distribution in the 25-150 keV range and is 29%, which 
indicates compatibility of our result with the background 
expectation. The inclusion of the neutron background 
has a negligible effect on the low-energy range due to the 
dominant surface-event background. 

We used the optimum interval method 2J| to com- 
pute the 90% C.L. upper limit on the spin- independent 
cross section as a function of WIMP mass and splitting. 
All eleven WIMP candidates were included as possible 
signal, with no background subtraction. The differen- 
tial rates were calculated under standard halo assump- 
tions according to (25j . which gives an updated version 
of the standard formula from 23|, correctly taking the 
effect of the Earth's velocity on the escape-velocity cut- 
off into account. We assumed this escape velocity Wcsc to 
be 544km/s [sij, while the standard value of 220km/s 
was applied for the dispersion vq of the Maxwellian dark 
matter velocity distribution. Helm form factors and a 
three-dimensional parametrization of the Earth velocity 
We were used following p3| . 

Regions allowed by DAM A/LIBRA at two different 
C.L.s (90,99.9%) were computed based on the published 
modulated spectrum in 6| from an exposure of 1.1 7 ton- 
years. As in [7}, we followed the goodness-of-fit tech- 
nique advocated in (25j to investigate the compatibility 
between the results from DAM A/LIBRA and CDMS. 
Quenching factors of 0.30 and 0.09 were applied for Na 
and I nuclei in the DAMA/LIBRA setup, respectively 
32 1 . Possible channeling effects [s^ were not included 
in this study since they do not have a significant impact 
on the results from an iDM analysis where a signal is 
expected at tens of kcV recoil energy 34 1 . 

Selected results from these computations are shown 
in Fig. |5] in the cross-section versus WIMP-mass plane 
for two chosen WIMP-mass splittings. The left plot 
shows the standard case with S = OkeV, equivalent 
to assuming elastic scattering, while 6 = 120 keV is 
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FIG. 8. (color online). 90% C.L. upper limits on the scalar WIMP-nucleon cross section for WIMP-mass splittings of keV (left) 
and 120 keV (right) from this analysis (red/dashed) and from our previous analysis (black/solid) The red/dotted line in the 
right plot indicates the expected sensitivity for this analysis based on our estimate of the total background. The colored regions 
represent DAMA/LIBRA allowed regions at two different C.L.s (90,99.9%) calculated following a ^ goodness-of-fit technique 
[25| . The cross ( x ) marks the parameter-space point which yields the minimum in the shown cross-section versus WIMP-mass 
plane given the WIMP-mass splitting. 



used for the right plot, a value vifhich is not experimen- 
tally excluded by our previous analysis. Apart from the 
DAMA/LIBRA allowed regions, and constraints emerg- 
ing from the analysis presented in this paper, the plots 
also contain cross-section limits from our previous analy- 
sis of the 10-100 keV energy range [7]. Constraints from 
the new analysis are less stringent. This was anticipated 
for the elastic scattering case and low WIMP-mass split- 
tings in general, since more surface-background events 
were expected at low energies due to the looser timing 
cut defined in the 25-150 keV energy range. The limits 
are slightly weaker at (5 = 120 keV, due to the occurrence 
of the three candidates above 25keV, where the rate is 
expected to peak for higher WIMP-mass splittings. The 
eight low-energy events have no effect on this part of the 
parameter space due to the utilization of the optimum 
interval method. WIMP masses above ^100 GeV/c^ are 
excluded for this mass splitting by the current and pre- 
vious analysis. 

Since the iDM parameter space is three-dimensional, 
consisting of the cross section, WIMP mass, and WIMP- 
mass splitting, we defined excluded regions by requiring 
the 90% C.L. upper limit on the cross section from CDMS 
to completely rule out the corresponding DAMA/LIBRA 
allowed cross sections (also at 90% C.L.) for given WIMP 
mass and WIMP-mass splitting. The results are shown 
in Fig. ini The only remaining allowed parameter space is 
within a narrow region at WIMP masses of ^100 GeV/c^ 
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FIG. 9. (color online). The blue/shaded regions represent 
WIMP masses and WIMP-mass splittings for which cross sec- 
tions exist that are compatible with the modulation spectrum 
observed by DAMA/LIBRA at 90% C.L. The hatched re- 
gions show constraints on these parameters from the analysis 
presented in this paper (red/dashed) and from our previous 
analysis (black/solid) [Tf. The black/dashed line represents 
the maximum reach of the CDMS II experiment. 
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FIG. 10. (color online). Constraints from CDMS on the iDM parameter space allowed by DAMA/LIBRA. Same as Fig. [9]but 
with different velocity-distribution parameters, as given in the plots. All other parameters remain unchanged. 



and WIMP-mass splittings between 85keV and 135 keV. 
In the case of the new analysis presented in this pa- 
per there is also a tiny area in the low-mass region 
which is not excluded. The black/dashed line repre- 
sents the maximum reach in the shown parameter space 
of an experiment using a Ge target like CDMS II. It 
is computed based entirely on kinematics by demand- 
ing Wmin = Vesc + "^E, and is therefore independent of 
the cross-section parameter. Even with higher exposure 
and increased sensitivity, CDMS II cannot rule out the 
entire DAMA/LIBRA allowed parameter space simply 



because (relative to an I nucleus) the Ge nucleus is too 
light. This is the main reason why the constraints from 
ZEPLIN-III [3^, which employs a Xe target, are more 
stringent. Nevertheless, the results from CDMS II cur- 
rently have competitive sensitivity compared to the con- 
straints from other Xe based experiments 3^ 3 



The iDM scenario with a nonzero b is particularly sen- 
sitive to the high- velocity tail of the dark matter velocity 
distribution due to the increased minimal velocity (see 
Eq. ([T])). Therefore, it exhibits a strong dependence on 
the velocity-distribution parameters vq (the dispersion) 



and Vesc (the galactic escape velocity) [3g]. To examine 
these dependencies, in Fig. [10] we show plots similar to 
Fig. [5] but with different values of Vcsc and Vq- The top 
plots explore the v^sc 90% C.L. lower and upper limits 
found in 3l| (498km/s and 608km/s), while all other 
parameters remain unchanged relative to Fig. [3] In the 
bottom plots we varied only vq, assigning a (convenient) 
lower value of 200 km/s for the left plot and a higher value 
of 254 km/s (the preferred value according to [39|) for the 
right plot. The capability of CDMS to constrain an iDM 
interpretation of the DAMA/LIBRA results is relatively 
independent of the actual velocity-distribution parame- 
ters. However, the shape and location of the parameter- 
space region, which is still allowed by CDMS, has a con- 
siderable dependence on Wcsc and vq, as expected. Non- 
Maxwellian velocity distributions as discussed in 
are beyond the scope of this study. 

In this paper we presented the first CDMS analysis 
which includes recoil energies up to 150 keV. The entire 
five-tower data set was used in a combined analysis. Be- 
cause of the occurrence of the three candidate events be- 
tween 25 keV and 150 keV the constraints on the iDM pa- 
rameter space are slightly weaker than from our previous 
analysis for which no events were observed at intermedi- 
ate energies where the rate is expected to peak. The only 
remaining parameter space allowed by CDMS data is 
within a narrow region at WIMP masses of ^100 GeV/c^ 
and WIMP-mass splittings between 85keV and 135 keV, 
assuming standard values for the WIMP-velocity distri- 
bution parameters. Varying the values of these parame- 
ters changes this region considerably but has only a mi- 
nor effect on the capability of CDMS to constrain an 
iDM interpretation of the DAMA/LIBRA results. Fi- 
nally, though this analysis was performed with regard to 
the iDM scenario, the expansion of the analysis range to 
150 keV could be useful to test other models predicting 
a signal at tens of keV recoil energy. 
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